\(\int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx\) [488]

Optimal result

Integrand size = 26, antiderivative size = 435 \[ \int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\frac {f x}{2 a d}-\frac {(e+f x)^2}{2 a f}-\frac {b^2 (e+f x)^2}{2 a^3 f}+\frac {\left (a^2+b^2\right ) (e+f x)^2}{2 a^3 f}+\frac {b f \text {arctanh}(\cosh (c+d x))}{a^2 d^2}-\frac {f \coth (c+d x)}{2 a d^2}-\frac {(e+f x) \coth ^2(c+d x)}{2 a d}+\frac {b (e+f x) \text {csch}(c+d x)}{a^2 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{a^3 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{a^3 d}+\frac {(e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a d}+\frac {b^2 (e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a^3 d}-\frac {\left (a^2+b^2\right ) f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{a^3 d^2}-\frac {\left (a^2+b^2\right ) f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{a^3 d^2}+\frac {f \operatorname {PolyLog}\left (2,e^{2 (c+d x)}\right )}{2 a d^2}+\frac {b^2 f \operatorname {PolyLog}\left (2,e^{2 (c+d x)}\right )}{2 a^3 d^2} \]

[Out]

1/2*f*x/a/d-1/2*(f*x+e)^2/a/f-1/2*b^2*(f*x+e)^2/a^3/f+1/2*(a^2+b^2)*(f*x+e)^2/a^3/f+b*f*arctanh(cosh(d*x+c))/a
^2/d^2-1/2*f*coth(d*x+c)/a/d^2-1/2*(f*x+e)*coth(d*x+c)^2/a/d+b*(f*x+e)*csch(d*x+c)/a^2/d+(f*x+e)*ln(1-exp(2*d*
x+2*c))/a/d+b^2*(f*x+e)*ln(1-exp(2*d*x+2*c))/a^3/d-(a^2+b^2)*(f*x+e)*ln(1+b*exp(d*x+c)/(a-(a^2+b^2)^(1/2)))/a^
3/d-(a^2+b^2)*(f*x+e)*ln(1+b*exp(d*x+c)/(a+(a^2+b^2)^(1/2)))/a^3/d+1/2*f*polylog(2,exp(2*d*x+2*c))/a/d^2+1/2*b
^2*f*polylog(2,exp(2*d*x+2*c))/a^3/d^2-(a^2+b^2)*f*polylog(2,-b*exp(d*x+c)/(a-(a^2+b^2)^(1/2)))/a^3/d^2-(a^2+b
^2)*f*polylog(2,-b*exp(d*x+c)/(a+(a^2+b^2)^(1/2)))/a^3/d^2

Rubi [A] (verified)

Time = 0.66 (sec) , antiderivative size = 435, normalized size of antiderivative = 1.00, number of steps used = 36, number of rules used = 18, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.692, Rules used = {5688, 3801, 3554, 8, 3797, 2221, 2317, 2438, 5704, 5558, 3377, 2718, 5560, 3855, 5554, 2715, 5684, 5680} \[ \int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\frac {b^2 f \operatorname {PolyLog}\left (2,e^{2 (c+d x)}\right )}{2 a^3 d^2}+\frac {b^2 (e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a^3 d}-\frac {b^2 (e+f x)^2}{2 a^3 f}+\frac {b f \text {arctanh}(\cosh (c+d x))}{a^2 d^2}+\frac {b (e+f x) \text {csch}(c+d x)}{a^2 d}-\frac {f \left (a^2+b^2\right ) \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{a^3 d^2}-\frac {f \left (a^2+b^2\right ) \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{a^3 d^2}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}+1\right )}{a^3 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (\frac {b e^{c+d x}}{\sqrt {a^2+b^2}+a}+1\right )}{a^3 d}+\frac {\left (a^2+b^2\right ) (e+f x)^2}{2 a^3 f}+\frac {f \operatorname {PolyLog}\left (2,e^{2 (c+d x)}\right )}{2 a d^2}-\frac {f \coth (c+d x)}{2 a d^2}+\frac {(e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a d}-\frac {(e+f x) \coth ^2(c+d x)}{2 a d}+\frac {f x}{2 a d}-\frac {(e+f x)^2}{2 a f} \]

[In]

Int[((e + f*x)*Coth[c + d*x]^3)/(a + b*Sinh[c + d*x]),x]

[Out]

(f*x)/(2*a*d) - (e + f*x)^2/(2*a*f) - (b^2*(e + f*x)^2)/(2*a^3*f) + ((a^2 + b^2)*(e + f*x)^2)/(2*a^3*f) + (b*f
*ArcTanh[Cosh[c + d*x]])/(a^2*d^2) - (f*Coth[c + d*x])/(2*a*d^2) - ((e + f*x)*Coth[c + d*x]^2)/(2*a*d) + (b*(e
+ f*x)*Csch[c + d*x])/(a^2*d) - ((a^2 + b^2)*(e + f*x)*Log[1 + (b*E^(c + d*x))/(a - Sqrt[a^2 + b^2])])/(a^3*d
) - ((a^2 + b^2)*(e + f*x)*Log[1 + (b*E^(c + d*x))/(a + Sqrt[a^2 + b^2])])/(a^3*d) + ((e + f*x)*Log[1 - E^(2*(
c + d*x))])/(a*d) + (b^2*(e + f*x)*Log[1 - E^(2*(c + d*x))])/(a^3*d) - ((a^2 + b^2)*f*PolyLog[2, -((b*E^(c + d
*x))/(a - Sqrt[a^2 + b^2]))])/(a^3*d^2) - ((a^2 + b^2)*f*PolyLog[2, -((b*E^(c + d*x))/(a + Sqrt[a^2 + b^2]))])
/(a^3*d^2) + (f*PolyLog[2, E^(2*(c + d*x))])/(2*a*d^2) + (b^2*f*PolyLog[2, E^(2*(c + d*x))])/(2*a^3*d^2)

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 2221

Int[(((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)*((c_.) + (d_.)*(x_))^(m_.))/((a_) + (b_.)*((F_)^((g_.)*((e_.) +
(f_.)*(x_))))^(n_.)), x_Symbol] :> Simp[((c + d*x)^m/(b*f*g*n*Log[F]))*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x]
- Dist[d*(m/(b*f*g*n*Log[F])), Int[(c + d*x)^(m - 1)*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x], x] /; FreeQ[{F,
a, b, c, d, e, f, g, n}, x] && IGtQ[m, 0]

Rule 2317

Int[Log[(a_) + (b_.)*((F_)^((e_.)*((c_.) + (d_.)*(x_))))^(n_.)], x_Symbol] :> Dist[1/(d*e*n*Log[F]), Subst[Int
[Log[a + b*x]/x, x], x, (F^(e*(c + d*x)))^n], x] /; FreeQ[{F, a, b, c, d, e, n}, x] && GtQ[a, 0]

Rule 2438

Int[Log[(c_.)*((d_) + (e_.)*(x_)^(n_.))]/(x_), x_Symbol] :> Simp[-PolyLog[2, (-c)*e*x^n]/n, x] /; FreeQ[{c, d,
e, n}, x] && EqQ[c*d, 1]

Rule 2715

Int[((b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[(-b)*Cos[c + d*x]*((b*Sin[c + d*x])^(n - 1)/(d*n))
, x] + Dist[b^2*((n - 1)/n), Int[(b*Sin[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] && Integ
erQ[2*n]

Rule 2718

Int[sin[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-Cos[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 3377

Int[((c_.) + (d_.)*(x_))^(m_.)*sin[(e_.) + (f_.)*(x_)], x_Symbol] :> Simp[(-(c + d*x)^m)*(Cos[e + f*x]/f), x]
+ Dist[d*(m/f), Int[(c + d*x)^(m - 1)*Cos[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && GtQ[m, 0]

Rule 3554

Int[((b_.)*tan[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[b*((b*Tan[c + d*x])^(n - 1)/(d*(n - 1))), x] - Dis
t[b^2, Int[(b*Tan[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1]

Rule 3797

Int[((c_.) + (d_.)*(x_))^(m_.)*tan[(e_.) + Pi*(k_.) + (Complex[0, fz_])*(f_.)*(x_)], x_Symbol] :> Simp[(-I)*((
c + d*x)^(m + 1)/(d*(m + 1))), x] + Dist[2*I, Int[((c + d*x)^m*(E^(2*((-I)*e + f*fz*x))/(1 + E^(2*((-I)*e + f*
fz*x))/E^(2*I*k*Pi))))/E^(2*I*k*Pi), x], x] /; FreeQ[{c, d, e, f, fz}, x] && IntegerQ[4*k] && IGtQ[m, 0]

Rule 3801

Int[((c_.) + (d_.)*(x_))^(m_.)*((b_.)*tan[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[b*(c + d*x)^m*((b*Tan[e
+ f*x])^(n - 1)/(f*(n - 1))), x] + (-Dist[b*d*(m/(f*(n - 1))), Int[(c + d*x)^(m - 1)*(b*Tan[e + f*x])^(n - 1)
, x], x] - Dist[b^2, Int[(c + d*x)^m*(b*Tan[e + f*x])^(n - 2), x], x]) /; FreeQ[{b, c, d, e, f}, x] && GtQ[n,
1] && GtQ[m, 0]

Rule 3855

Int[csc[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-ArcTanh[Cos[c + d*x]]/d, x] /; FreeQ[{c, d}, x]

Rule 5554

Int[Cosh[(a_.) + (b_.)*(x_)]*((c_.) + (d_.)*(x_))^(m_.)*Sinh[(a_.) + (b_.)*(x_)]^(n_.), x_Symbol] :> Simp[(c +
d*x)^m*(Sinh[a + b*x]^(n + 1)/(b*(n + 1))), x] - Dist[d*(m/(b*(n + 1))), Int[(c + d*x)^(m - 1)*Sinh[a + b*x]^
(n + 1), x], x] /; FreeQ[{a, b, c, d, n}, x] && IGtQ[m, 0] && NeQ[n, -1]

Rule 5558

Int[Cosh[(a_.) + (b_.)*(x_)]^(n_.)*Coth[(a_.) + (b_.)*(x_)]^(p_.)*((c_.) + (d_.)*(x_))^(m_.), x_Symbol] :> Int
[(c + d*x)^m*Cosh[a + b*x]^n*Coth[a + b*x]^(p - 2), x] + Int[(c + d*x)^m*Cosh[a + b*x]^(n - 2)*Coth[a + b*x]^p
, x] /; FreeQ[{a, b, c, d, m}, x] && IGtQ[n, 0] && IGtQ[p, 0]

Rule 5560

Int[Coth[(a_.) + (b_.)*(x_)]^(p_.)*Csch[(a_.) + (b_.)*(x_)]^(n_.)*((c_.) + (d_.)*(x_))^(m_.), x_Symbol] :> Sim
p[(-(c + d*x)^m)*(Csch[a + b*x]^n/(b*n)), x] + Dist[d*(m/(b*n)), Int[(c + d*x)^(m - 1)*Csch[a + b*x]^n, x], x]
/; FreeQ[{a, b, c, d, n}, x] && EqQ[p, 1] && GtQ[m, 0]

Rule 5680

Int[(Cosh[(c_.) + (d_.)*(x_)]*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_.)*Sinh[(c_.) + (d_.)*(x_)]), x_Symbol] :
> Simp[-(e + f*x)^(m + 1)/(b*f*(m + 1)), x] + (Int[(e + f*x)^m*(E^(c + d*x)/(a - Rt[a^2 + b^2, 2] + b*E^(c + d
*x))), x] + Int[(e + f*x)^m*(E^(c + d*x)/(a + Rt[a^2 + b^2, 2] + b*E^(c + d*x))), x]) /; FreeQ[{a, b, c, d, e,
f}, x] && IGtQ[m, 0] && NeQ[a^2 + b^2, 0]

Rule 5684

Int[(Cosh[(c_.) + (d_.)*(x_)]^(n_)*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_.)*Sinh[(c_.) + (d_.)*(x_)]), x_Symb
ol] :> Dist[-a/b^2, Int[(e + f*x)^m*Cosh[c + d*x]^(n - 2), x], x] + (Dist[1/b, Int[(e + f*x)^m*Cosh[c + d*x]^(
n - 2)*Sinh[c + d*x], x], x] + Dist[(a^2 + b^2)/b^2, Int[(e + f*x)^m*(Cosh[c + d*x]^(n - 2)/(a + b*Sinh[c + d*
x])), x], x]) /; FreeQ[{a, b, c, d, e, f}, x] && IGtQ[n, 1] && NeQ[a^2 + b^2, 0] && IGtQ[m, 0]

Rule 5688

Int[(Coth[(c_.) + (d_.)*(x_)]^(n_.)*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_.)*Sinh[(c_.) + (d_.)*(x_)]), x_Sym
bol] :> Dist[1/a, Int[(e + f*x)^m*Coth[c + d*x]^n, x], x] - Dist[b/a, Int[(e + f*x)^m*Cosh[c + d*x]*(Coth[c +
d*x]^(n - 1)/(a + b*Sinh[c + d*x])), x], x] /; FreeQ[{a, b, c, d, e, f}, x] && IGtQ[m, 0] && IGtQ[n, 0]

Rule 5704

Int[(Cosh[(c_.) + (d_.)*(x_)]^(p_.)*Coth[(c_.) + (d_.)*(x_)]^(n_.)*((e_.) + (f_.)*(x_))^(m_.))/((a_) + (b_.)*S
inh[(c_.) + (d_.)*(x_)]), x_Symbol] :> Dist[1/a, Int[(e + f*x)^m*Cosh[c + d*x]^p*Coth[c + d*x]^n, x], x] - Dis
t[b/a, Int[(e + f*x)^m*Cosh[c + d*x]^(p + 1)*(Coth[c + d*x]^(n - 1)/(a + b*Sinh[c + d*x])), x], x] /; FreeQ[{a
, b, c, d, e, f}, x] && IGtQ[m, 0] && IGtQ[n, 0] && IGtQ[p, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\int (e+f x) \coth ^3(c+d x) \, dx}{a}-\frac {b \int \frac {(e+f x) \cosh (c+d x) \coth ^2(c+d x)}{a+b \sinh (c+d x)} \, dx}{a} \\ & = -\frac {(e+f x) \coth ^2(c+d x)}{2 a d}+\frac {\int (e+f x) \coth (c+d x) \, dx}{a}-\frac {b \int (e+f x) \cosh (c+d x) \coth ^2(c+d x) \, dx}{a^2}+\frac {b^2 \int \frac {(e+f x) \cosh ^2(c+d x) \coth (c+d x)}{a+b \sinh (c+d x)} \, dx}{a^2}+\frac {f \int \coth ^2(c+d x) \, dx}{2 a d} \\ & = -\frac {(e+f x)^2}{2 a f}-\frac {f \coth (c+d x)}{2 a d^2}-\frac {(e+f x) \coth ^2(c+d x)}{2 a d}-\frac {2 \int \frac {e^{2 (c+d x)} (e+f x)}{1-e^{2 (c+d x)}} \, dx}{a}-\frac {b \int (e+f x) \cosh (c+d x) \, dx}{a^2}-\frac {b \int (e+f x) \coth (c+d x) \text {csch}(c+d x) \, dx}{a^2}+\frac {b^2 \int (e+f x) \cosh ^2(c+d x) \coth (c+d x) \, dx}{a^3}-\frac {b^3 \int \frac {(e+f x) \cosh ^3(c+d x)}{a+b \sinh (c+d x)} \, dx}{a^3}+\frac {f \int 1 \, dx}{2 a d} \\ & = \frac {f x}{2 a d}-\frac {(e+f x)^2}{2 a f}-\frac {f \coth (c+d x)}{2 a d^2}-\frac {(e+f x) \coth ^2(c+d x)}{2 a d}+\frac {b (e+f x) \text {csch}(c+d x)}{a^2 d}+\frac {(e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a d}-\frac {b (e+f x) \sinh (c+d x)}{a^2 d}+\frac {b \int (e+f x) \cosh (c+d x) \, dx}{a^2}+\frac {b^2 \int (e+f x) \coth (c+d x) \, dx}{a^3}-\frac {\left (b \left (a^2+b^2\right )\right ) \int \frac {(e+f x) \cosh (c+d x)}{a+b \sinh (c+d x)} \, dx}{a^3}-\frac {f \int \log \left (1-e^{2 (c+d x)}\right ) \, dx}{a d}-\frac {(b f) \int \text {csch}(c+d x) \, dx}{a^2 d}+\frac {(b f) \int \sinh (c+d x) \, dx}{a^2 d} \\ & = \frac {f x}{2 a d}-\frac {(e+f x)^2}{2 a f}-\frac {b^2 (e+f x)^2}{2 a^3 f}+\frac {\left (a^2+b^2\right ) (e+f x)^2}{2 a^3 f}+\frac {b f \text {arctanh}(\cosh (c+d x))}{a^2 d^2}+\frac {b f \cosh (c+d x)}{a^2 d^2}-\frac {f \coth (c+d x)}{2 a d^2}-\frac {(e+f x) \coth ^2(c+d x)}{2 a d}+\frac {b (e+f x) \text {csch}(c+d x)}{a^2 d}+\frac {(e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a d}-\frac {\left (2 b^2\right ) \int \frac {e^{2 (c+d x)} (e+f x)}{1-e^{2 (c+d x)}} \, dx}{a^3}-\frac {\left (b \left (a^2+b^2\right )\right ) \int \frac {e^{c+d x} (e+f x)}{a-\sqrt {a^2+b^2}+b e^{c+d x}} \, dx}{a^3}-\frac {\left (b \left (a^2+b^2\right )\right ) \int \frac {e^{c+d x} (e+f x)}{a+\sqrt {a^2+b^2}+b e^{c+d x}} \, dx}{a^3}-\frac {f \text {Subst}\left (\int \frac {\log (1-x)}{x} \, dx,x,e^{2 (c+d x)}\right )}{2 a d^2}-\frac {(b f) \int \sinh (c+d x) \, dx}{a^2 d} \\ & = \frac {f x}{2 a d}-\frac {(e+f x)^2}{2 a f}-\frac {b^2 (e+f x)^2}{2 a^3 f}+\frac {\left (a^2+b^2\right ) (e+f x)^2}{2 a^3 f}+\frac {b f \text {arctanh}(\cosh (c+d x))}{a^2 d^2}-\frac {f \coth (c+d x)}{2 a d^2}-\frac {(e+f x) \coth ^2(c+d x)}{2 a d}+\frac {b (e+f x) \text {csch}(c+d x)}{a^2 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{a^3 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{a^3 d}+\frac {(e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a d}+\frac {b^2 (e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a^3 d}+\frac {f \operatorname {PolyLog}\left (2,e^{2 (c+d x)}\right )}{2 a d^2}-\frac {\left (b^2 f\right ) \int \log \left (1-e^{2 (c+d x)}\right ) \, dx}{a^3 d}+\frac {\left (\left (a^2+b^2\right ) f\right ) \int \log \left (1+\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right ) \, dx}{a^3 d}+\frac {\left (\left (a^2+b^2\right ) f\right ) \int \log \left (1+\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right ) \, dx}{a^3 d} \\ & = \frac {f x}{2 a d}-\frac {(e+f x)^2}{2 a f}-\frac {b^2 (e+f x)^2}{2 a^3 f}+\frac {\left (a^2+b^2\right ) (e+f x)^2}{2 a^3 f}+\frac {b f \text {arctanh}(\cosh (c+d x))}{a^2 d^2}-\frac {f \coth (c+d x)}{2 a d^2}-\frac {(e+f x) \coth ^2(c+d x)}{2 a d}+\frac {b (e+f x) \text {csch}(c+d x)}{a^2 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{a^3 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{a^3 d}+\frac {(e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a d}+\frac {b^2 (e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a^3 d}+\frac {f \operatorname {PolyLog}\left (2,e^{2 (c+d x)}\right )}{2 a d^2}-\frac {\left (b^2 f\right ) \text {Subst}\left (\int \frac {\log (1-x)}{x} \, dx,x,e^{2 (c+d x)}\right )}{2 a^3 d^2}+\frac {\left (\left (a^2+b^2\right ) f\right ) \text {Subst}\left (\int \frac {\log \left (1+\frac {b x}{a-\sqrt {a^2+b^2}}\right )}{x} \, dx,x,e^{c+d x}\right )}{a^3 d^2}+\frac {\left (\left (a^2+b^2\right ) f\right ) \text {Subst}\left (\int \frac {\log \left (1+\frac {b x}{a+\sqrt {a^2+b^2}}\right )}{x} \, dx,x,e^{c+d x}\right )}{a^3 d^2} \\ & = \frac {f x}{2 a d}-\frac {(e+f x)^2}{2 a f}-\frac {b^2 (e+f x)^2}{2 a^3 f}+\frac {\left (a^2+b^2\right ) (e+f x)^2}{2 a^3 f}+\frac {b f \text {arctanh}(\cosh (c+d x))}{a^2 d^2}-\frac {f \coth (c+d x)}{2 a d^2}-\frac {(e+f x) \coth ^2(c+d x)}{2 a d}+\frac {b (e+f x) \text {csch}(c+d x)}{a^2 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{a^3 d}-\frac {\left (a^2+b^2\right ) (e+f x) \log \left (1+\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{a^3 d}+\frac {(e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a d}+\frac {b^2 (e+f x) \log \left (1-e^{2 (c+d x)}\right )}{a^3 d}-\frac {\left (a^2+b^2\right ) f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )}{a^3 d^2}-\frac {\left (a^2+b^2\right ) f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )}{a^3 d^2}+\frac {f \operatorname {PolyLog}\left (2,e^{2 (c+d x)}\right )}{2 a d^2}+\frac {b^2 f \operatorname {PolyLog}\left (2,e^{2 (c+d x)}\right )}{2 a^3 d^2} \\ \end{align*}

Mathematica [A] (warning: unable to verify)

Time = 8.63 (sec) , antiderivative size = 766, normalized size of antiderivative = 1.76 \[ \int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\frac {\left (2 b d e \cosh \left (\frac {1}{2} (c+d x)\right )-a f \cosh \left (\frac {1}{2} (c+d x)\right )-2 b c f \cosh \left (\frac {1}{2} (c+d x)\right )+2 b f (c+d x) \cosh \left (\frac {1}{2} (c+d x)\right )\right ) \text {csch}\left (\frac {1}{2} (c+d x)\right )}{4 a^2 d^2}+\frac {(-d e+c f-f (c+d x)) \text {csch}^2\left (\frac {1}{2} (c+d x)\right )}{8 a d^2}+\frac {\frac {\left (a^2+b^2\right ) (d e+d f x)^2}{2 f}+\left (-a b f+a^2 (d e+d f x)+b^2 (d e+d f x)\right ) \log \left (1-e^{-c-d x}\right )+\left (a b f+a^2 (d e+d f x)+b^2 (d e+d f x)\right ) \log \left (1+e^{-c-d x}\right )-\left (a^2+b^2\right ) f \operatorname {PolyLog}\left (2,-e^{-c-d x}\right )-\left (a^2+b^2\right ) f \operatorname {PolyLog}\left (2,e^{-c-d x}\right )}{a^3 d^2}-\frac {\left (a^2+b^2\right ) \left (-2 d e (c+d x)+2 c f (c+d x)-f (c+d x)^2+\frac {4 a \sqrt {a^2+b^2} d e \arctan \left (\frac {a+b e^{c+d x}}{\sqrt {-a^2-b^2}}\right )}{\sqrt {-\left (a^2+b^2\right )^2}}-\frac {4 a \sqrt {-\left (a^2+b^2\right )^2} d e \text {arctanh}\left (\frac {a+b e^{c+d x}}{\sqrt {a^2+b^2}}\right )}{\left (-a^2-b^2\right )^{3/2}}+2 f (c+d x) \log \left (1+\frac {b e^{c+d x}}{a-\sqrt {a^2+b^2}}\right )+2 f (c+d x) \log \left (1+\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )-2 c f \log \left (b-2 a e^{c+d x}-b e^{2 (c+d x)}\right )+2 d e \log \left (2 a e^{c+d x}+b \left (-1+e^{2 (c+d x)}\right )\right )+2 f \operatorname {PolyLog}\left (2,\frac {b e^{c+d x}}{-a+\sqrt {a^2+b^2}}\right )+2 f \operatorname {PolyLog}\left (2,-\frac {b e^{c+d x}}{a+\sqrt {a^2+b^2}}\right )\right )}{2 a^3 d^2}+\frac {(d e-c f+f (c+d x)) \text {sech}^2\left (\frac {1}{2} (c+d x)\right )}{8 a d^2}+\frac {\text {sech}\left (\frac {1}{2} (c+d x)\right ) \left (-2 b d e \sinh \left (\frac {1}{2} (c+d x)\right )-a f \sinh \left (\frac {1}{2} (c+d x)\right )+2 b c f \sinh \left (\frac {1}{2} (c+d x)\right )-2 b f (c+d x) \sinh \left (\frac {1}{2} (c+d x)\right )\right )}{4 a^2 d^2} \]

[In]

Integrate[((e + f*x)*Coth[c + d*x]^3)/(a + b*Sinh[c + d*x]),x]

[Out]

((2*b*d*e*Cosh[(c + d*x)/2] - a*f*Cosh[(c + d*x)/2] - 2*b*c*f*Cosh[(c + d*x)/2] + 2*b*f*(c + d*x)*Cosh[(c + d*
x)/2])*Csch[(c + d*x)/2])/(4*a^2*d^2) + ((-(d*e) + c*f - f*(c + d*x))*Csch[(c + d*x)/2]^2)/(8*a*d^2) + (((a^2
+ b^2)*(d*e + d*f*x)^2)/(2*f) + (-(a*b*f) + a^2*(d*e + d*f*x) + b^2*(d*e + d*f*x))*Log[1 - E^(-c - d*x)] + (a*
b*f + a^2*(d*e + d*f*x) + b^2*(d*e + d*f*x))*Log[1 + E^(-c - d*x)] - (a^2 + b^2)*f*PolyLog[2, -E^(-c - d*x)] -
(a^2 + b^2)*f*PolyLog[2, E^(-c - d*x)])/(a^3*d^2) - ((a^2 + b^2)*(-2*d*e*(c + d*x) + 2*c*f*(c + d*x) - f*(c +
d*x)^2 + (4*a*Sqrt[a^2 + b^2]*d*e*ArcTan[(a + b*E^(c + d*x))/Sqrt[-a^2 - b^2]])/Sqrt[-(a^2 + b^2)^2] - (4*a*S
qrt[-(a^2 + b^2)^2]*d*e*ArcTanh[(a + b*E^(c + d*x))/Sqrt[a^2 + b^2]])/(-a^2 - b^2)^(3/2) + 2*f*(c + d*x)*Log[1
+ (b*E^(c + d*x))/(a - Sqrt[a^2 + b^2])] + 2*f*(c + d*x)*Log[1 + (b*E^(c + d*x))/(a + Sqrt[a^2 + b^2])] - 2*c
*f*Log[b - 2*a*E^(c + d*x) - b*E^(2*(c + d*x))] + 2*d*e*Log[2*a*E^(c + d*x) + b*(-1 + E^(2*(c + d*x)))] + 2*f*
PolyLog[2, (b*E^(c + d*x))/(-a + Sqrt[a^2 + b^2])] + 2*f*PolyLog[2, -((b*E^(c + d*x))/(a + Sqrt[a^2 + b^2]))])
)/(2*a^3*d^2) + ((d*e - c*f + f*(c + d*x))*Sech[(c + d*x)/2]^2)/(8*a*d^2) + (Sech[(c + d*x)/2]*(-2*b*d*e*Sinh[
(c + d*x)/2] - a*f*Sinh[(c + d*x)/2] + 2*b*c*f*Sinh[(c + d*x)/2] - 2*b*f*(c + d*x)*Sinh[(c + d*x)/2]))/(4*a^2*
d^2)

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(1097\) vs. \(2(407)=814\).

Time = 1.98 (sec) , antiderivative size = 1098, normalized size of antiderivative = 2.52

method	result	size
risch	\(\text {Expression too large to display}\)	\(1098\)

[In]

int((f*x+e)*coth(d*x+c)^3/(a+b*sinh(d*x+c)),x,method=_RETURNVERBOSE)

[Out]

-(-2*b*d*f*x*exp(3*d*x+3*c)+2*a*d*f*x*exp(2*d*x+2*c)-2*b*d*e*exp(3*d*x+3*c)+2*a*d*e*exp(2*d*x+2*c)+2*b*d*f*x*e
xp(d*x+c)+a*f*exp(2*d*x+2*c)+2*b*d*e*exp(d*x+c)-a*f)/a^2/d^2/(exp(2*d*x+2*c)-1)^2-1/d*f/a*ln((b*exp(d*x+c)+(a^
2+b^2)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))*x-1/d*f/a*ln((-b*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))*x+1/
d*f/a*ln(exp(d*x+c)+1)*x-1/d^2*f/a*ln((b*exp(d*x+c)+(a^2+b^2)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))*c-1/d^2*f/a*ln((-b
*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))*c-1/d^2*c*f/a*ln(exp(d*x+c)-1)+1/d^2*c*f/a*ln(b*exp(2*d*x
+2*c)+2*a*exp(d*x+c)-b)-1/d^2*b^2/a^3*f*dilog(exp(d*x+c))+1/d^2*b^2/a^3*f*dilog(exp(d*x+c)+1)+1/d*b^2/a^3*e*ln
(exp(d*x+c)-1)+1/d*b^2/a^3*e*ln(exp(d*x+c)+1)-1/d*b^2/a^3*e*ln(b*exp(2*d*x+2*c)+2*a*exp(d*x+c)-b)-1/d^2*b^2/a^
3*f*dilog((-b*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))-1/d^2*b^2/a^3*f*dilog((b*exp(d*x+c)+(a^2+b^2
)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))-1/d^2*b/a^2*f*ln(exp(d*x+c)-1)+1/d^2*b/a^2*f*ln(exp(d*x+c)+1)-1/d^2*f/a*dilog(
(-b*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))-1/d^2*f/a*dilog((b*exp(d*x+c)+(a^2+b^2)^(1/2)+a)/(a+(a
^2+b^2)^(1/2)))-1/d^2*f/a*dilog(exp(d*x+c))+1/d*e/a*ln(exp(d*x+c)-1)-1/d^2*b^2/a^3*f*ln((b*exp(d*x+c)+(a^2+b^2
)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))*c-1/d^2*b^2/a^3*f*ln((-b*exp(d*x+c)+(a^2+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))*c
-1/d*b^2/a^3*f*ln((b*exp(d*x+c)+(a^2+b^2)^(1/2)+a)/(a+(a^2+b^2)^(1/2)))*x-1/d*b^2/a^3*f*ln((-b*exp(d*x+c)+(a^2
+b^2)^(1/2)-a)/(-a+(a^2+b^2)^(1/2)))*x+1/d*e/a*ln(exp(d*x+c)+1)-1/d*e/a*ln(b*exp(2*d*x+2*c)+2*a*exp(d*x+c)-b)+
1/d^2*f/a*dilog(exp(d*x+c)+1)+1/d*b^2/a^3*f*ln(exp(d*x+c)+1)*x-1/d^2*b^2/a^3*c*f*ln(exp(d*x+c)-1)+1/d^2*b^2/a^
3*c*f*ln(b*exp(2*d*x+2*c)+2*a*exp(d*x+c)-b)

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 3547 vs. \(2 (403) = 806\).

Time = 0.31 (sec) , antiderivative size = 3547, normalized size of antiderivative = 8.15 \[ \int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\text {Too large to display} \]

[In]

integrate((f*x+e)*coth(d*x+c)^3/(a+b*sinh(d*x+c)),x, algorithm="fricas")

[Out]

(2*(a*b*d*f*x + a*b*d*e)*cosh(d*x + c)^3 + 2*(a*b*d*f*x + a*b*d*e)*sinh(d*x + c)^3 + a^2*f - (2*a^2*d*f*x + 2*
a^2*d*e + a^2*f)*cosh(d*x + c)^2 - (2*a^2*d*f*x + 2*a^2*d*e + a^2*f - 6*(a*b*d*f*x + a*b*d*e)*cosh(d*x + c))*s
inh(d*x + c)^2 - 2*(a*b*d*f*x + a*b*d*e)*cosh(d*x + c) - ((a^2 + b^2)*f*cosh(d*x + c)^4 + 4*(a^2 + b^2)*f*cosh
(d*x + c)*sinh(d*x + c)^3 + (a^2 + b^2)*f*sinh(d*x + c)^4 - 2*(a^2 + b^2)*f*cosh(d*x + c)^2 + 2*(3*(a^2 + b^2)
*f*cosh(d*x + c)^2 - (a^2 + b^2)*f)*sinh(d*x + c)^2 + (a^2 + b^2)*f + 4*((a^2 + b^2)*f*cosh(d*x + c)^3 - (a^2
+ b^2)*f*cosh(d*x + c))*sinh(d*x + c))*dilog((a*cosh(d*x + c) + a*sinh(d*x + c) + (b*cosh(d*x + c) + b*sinh(d*
x + c))*sqrt((a^2 + b^2)/b^2) - b)/b + 1) - ((a^2 + b^2)*f*cosh(d*x + c)^4 + 4*(a^2 + b^2)*f*cosh(d*x + c)*sin
h(d*x + c)^3 + (a^2 + b^2)*f*sinh(d*x + c)^4 - 2*(a^2 + b^2)*f*cosh(d*x + c)^2 + 2*(3*(a^2 + b^2)*f*cosh(d*x +
c)^2 - (a^2 + b^2)*f)*sinh(d*x + c)^2 + (a^2 + b^2)*f + 4*((a^2 + b^2)*f*cosh(d*x + c)^3 - (a^2 + b^2)*f*cosh
(d*x + c))*sinh(d*x + c))*dilog((a*cosh(d*x + c) + a*sinh(d*x + c) - (b*cosh(d*x + c) + b*sinh(d*x + c))*sqrt(
(a^2 + b^2)/b^2) - b)/b + 1) + ((a^2 + b^2)*f*cosh(d*x + c)^4 + 4*(a^2 + b^2)*f*cosh(d*x + c)*sinh(d*x + c)^3
+ (a^2 + b^2)*f*sinh(d*x + c)^4 - 2*(a^2 + b^2)*f*cosh(d*x + c)^2 + 2*(3*(a^2 + b^2)*f*cosh(d*x + c)^2 - (a^2
+ b^2)*f)*sinh(d*x + c)^2 + (a^2 + b^2)*f + 4*((a^2 + b^2)*f*cosh(d*x + c)^3 - (a^2 + b^2)*f*cosh(d*x + c))*si
nh(d*x + c))*dilog(cosh(d*x + c) + sinh(d*x + c)) + ((a^2 + b^2)*f*cosh(d*x + c)^4 + 4*(a^2 + b^2)*f*cosh(d*x
+ c)*sinh(d*x + c)^3 + (a^2 + b^2)*f*sinh(d*x + c)^4 - 2*(a^2 + b^2)*f*cosh(d*x + c)^2 + 2*(3*(a^2 + b^2)*f*co
sh(d*x + c)^2 - (a^2 + b^2)*f)*sinh(d*x + c)^2 + (a^2 + b^2)*f + 4*((a^2 + b^2)*f*cosh(d*x + c)^3 - (a^2 + b^2
)*f*cosh(d*x + c))*sinh(d*x + c))*dilog(-cosh(d*x + c) - sinh(d*x + c)) - (((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)
*cosh(d*x + c)^4 + 4*((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)*sinh(d*x + c)^3 + ((a^2 + b^2)*d*e - (a
^2 + b^2)*c*f)*sinh(d*x + c)^4 + (a^2 + b^2)*d*e - (a^2 + b^2)*c*f - 2*((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cos
h(d*x + c)^2 - 2*((a^2 + b^2)*d*e - (a^2 + b^2)*c*f - 3*((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)^2)*s
inh(d*x + c)^2 + 4*(((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)^3 - ((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*
cosh(d*x + c))*sinh(d*x + c))*log(2*b*cosh(d*x + c) + 2*b*sinh(d*x + c) + 2*b*sqrt((a^2 + b^2)/b^2) + 2*a) - (
((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)^4 + 4*((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)*sinh
(d*x + c)^3 + ((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*sinh(d*x + c)^4 + (a^2 + b^2)*d*e - (a^2 + b^2)*c*f - 2*((a^
2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)^2 - 2*((a^2 + b^2)*d*e - (a^2 + b^2)*c*f - 3*((a^2 + b^2)*d*e -
(a^2 + b^2)*c*f)*cosh(d*x + c)^2)*sinh(d*x + c)^2 + 4*(((a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c)^3 - (
(a^2 + b^2)*d*e - (a^2 + b^2)*c*f)*cosh(d*x + c))*sinh(d*x + c))*log(2*b*cosh(d*x + c) + 2*b*sinh(d*x + c) - 2
*b*sqrt((a^2 + b^2)/b^2) + 2*a) - (((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^4 + 4*((a^2 + b^2)*d*f*
x + (a^2 + b^2)*c*f)*cosh(d*x + c)*sinh(d*x + c)^3 + ((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*sinh(d*x + c)^4 + (
a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f - 2*((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^2 - 2*((a^2 + b^2)*
d*f*x + (a^2 + b^2)*c*f - 3*((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^2)*sinh(d*x + c)^2 + 4*(((a^2
+ b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^3 - ((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c))*sinh(d*
x + c))*log(-(a*cosh(d*x + c) + a*sinh(d*x + c) + (b*cosh(d*x + c) + b*sinh(d*x + c))*sqrt((a^2 + b^2)/b^2) -
b)/b) - (((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^4 + 4*((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(
d*x + c)*sinh(d*x + c)^3 + ((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*sinh(d*x + c)^4 + (a^2 + b^2)*d*f*x + (a^2 +
b^2)*c*f - 2*((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^2 - 2*((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f -
3*((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^2)*sinh(d*x + c)^2 + 4*(((a^2 + b^2)*d*f*x + (a^2 + b^2)
*c*f)*cosh(d*x + c)^3 - ((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c))*sinh(d*x + c))*log(-(a*cosh(d*x +
c) + a*sinh(d*x + c) - (b*cosh(d*x + c) + b*sinh(d*x + c))*sqrt((a^2 + b^2)/b^2) - b)/b) + (((a^2 + b^2)*d*f*
x + (a^2 + b^2)*d*e + a*b*f)*cosh(d*x + c)^4 + 4*((a^2 + b^2)*d*f*x + (a^2 + b^2)*d*e + a*b*f)*cosh(d*x + c)*s
inh(d*x + c)^3 + ((a^2 + b^2)*d*f*x + (a^2 + b^2)*d*e + a*b*f)*sinh(d*x + c)^4 + (a^2 + b^2)*d*f*x + (a^2 + b^
2)*d*e + a*b*f - 2*((a^2 + b^2)*d*f*x + (a^2 + b^2)*d*e + a*b*f)*cosh(d*x + c)^2 - 2*((a^2 + b^2)*d*f*x + (a^2
+ b^2)*d*e + a*b*f - 3*((a^2 + b^2)*d*f*x + (a^2 + b^2)*d*e + a*b*f)*cosh(d*x + c)^2)*sinh(d*x + c)^2 + 4*(((
a^2 + b^2)*d*f*x + (a^2 + b^2)*d*e + a*b*f)*cosh(d*x + c)^3 - ((a^2 + b^2)*d*f*x + (a^2 + b^2)*d*e + a*b*f)*co
sh(d*x + c))*sinh(d*x + c))*log(cosh(d*x + c) + sinh(d*x + c) + 1) + (((a^2 + b^2)*d*e - (a*b + (a^2 + b^2)*c)
*f)*cosh(d*x + c)^4 + 4*((a^2 + b^2)*d*e - (a*b + (a^2 + b^2)*c)*f)*cosh(d*x + c)*sinh(d*x + c)^3 + ((a^2 + b^
2)*d*e - (a*b + (a^2 + b^2)*c)*f)*sinh(d*x + c)^4 + (a^2 + b^2)*d*e - 2*((a^2 + b^2)*d*e - (a*b + (a^2 + b^2)*
c)*f)*cosh(d*x + c)^2 - 2*((a^2 + b^2)*d*e - 3*((a^2 + b^2)*d*e - (a*b + (a^2 + b^2)*c)*f)*cosh(d*x + c)^2 - (
a*b + (a^2 + b^2)*c)*f)*sinh(d*x + c)^2 - (a*b + (a^2 + b^2)*c)*f + 4*(((a^2 + b^2)*d*e - (a*b + (a^2 + b^2)*c
)*f)*cosh(d*x + c)^3 - ((a^2 + b^2)*d*e - (a*b + (a^2 + b^2)*c)*f)*cosh(d*x + c))*sinh(d*x + c))*log(cosh(d*x
+ c) + sinh(d*x + c) - 1) + (((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^4 + 4*((a^2 + b^2)*d*f*x + (a
^2 + b^2)*c*f)*cosh(d*x + c)*sinh(d*x + c)^3 + ((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*sinh(d*x + c)^4 + (a^2 +
b^2)*d*f*x + (a^2 + b^2)*c*f - 2*((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^2 - 2*((a^2 + b^2)*d*f*x
+ (a^2 + b^2)*c*f - 3*((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^2)*sinh(d*x + c)^2 + 4*(((a^2 + b^2)
*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c)^3 - ((a^2 + b^2)*d*f*x + (a^2 + b^2)*c*f)*cosh(d*x + c))*sinh(d*x + c)
)*log(-cosh(d*x + c) - sinh(d*x + c) + 1) - 2*(a*b*d*f*x + a*b*d*e - 3*(a*b*d*f*x + a*b*d*e)*cosh(d*x + c)^2 +
(2*a^2*d*f*x + 2*a^2*d*e + a^2*f)*cosh(d*x + c))*sinh(d*x + c))/(a^3*d^2*cosh(d*x + c)^4 + 4*a^3*d^2*cosh(d*x
+ c)*sinh(d*x + c)^3 + a^3*d^2*sinh(d*x + c)^4 - 2*a^3*d^2*cosh(d*x + c)^2 + a^3*d^2 + 2*(3*a^3*d^2*cosh(d*x
+ c)^2 - a^3*d^2)*sinh(d*x + c)^2 + 4*(a^3*d^2*cosh(d*x + c)^3 - a^3*d^2*cosh(d*x + c))*sinh(d*x + c))

Sympy [F]

\[ \int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\int \frac {\left (e + f x\right ) \coth ^{3}{\left (c + d x \right )}}{a + b \sinh {\left (c + d x \right )}}\, dx \]

[In]

integrate((f*x+e)*coth(d*x+c)**3/(a+b*sinh(d*x+c)),x)

[Out]

Integral((e + f*x)*coth(c + d*x)**3/(a + b*sinh(c + d*x)), x)

Maxima [F]

\[ \int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\int { \frac {{\left (f x + e\right )} \coth \left (d x + c\right )^{3}}{b \sinh \left (d x + c\right ) + a} \,d x } \]

[In]

integrate((f*x+e)*coth(d*x+c)^3/(a+b*sinh(d*x+c)),x, algorithm="maxima")

[Out]

-(a^2*d*integrate(x/(a^3*d*e^(d*x + c) + a^3*d), x) + b^2*d*integrate(x/(a^3*d*e^(d*x + c) + a^3*d), x) - a^2*
d*integrate(x/(a^3*d*e^(d*x + c) - a^3*d), x) - b^2*d*integrate(x/(a^3*d*e^(d*x + c) - a^3*d), x) + a*b*((d*x
+ c)/(a^3*d^2) - log(e^(d*x + c) + 1)/(a^3*d^2)) - a*b*((d*x + c)/(a^3*d^2) - log(e^(d*x + c) - 1)/(a^3*d^2))
- (2*b*d*x*e^(3*d*x + 3*c) - 2*b*d*x*e^(d*x + c) - (2*a*d*x*e^(2*c) + a*e^(2*c))*e^(2*d*x) + a)/(a^2*d^2*e^(4*
d*x + 4*c) - 2*a^2*d^2*e^(2*d*x + 2*c) + a^2*d^2) - integrate(2*((a^3*e^c + a*b^2*e^c)*x*e^(d*x) - (a^2*b + b^
3)*x)/(a^3*b*e^(2*d*x + 2*c) + 2*a^4*e^(d*x + c) - a^3*b), x))*f - e*(2*(b*e^(-d*x - c) - a*e^(-2*d*x - 2*c) -
b*e^(-3*d*x - 3*c))/((2*a^2*e^(-2*d*x - 2*c) - a^2*e^(-4*d*x - 4*c) - a^2)*d) + (a^2 + b^2)*log(-2*a*e^(-d*x
- c) + b*e^(-2*d*x - 2*c) - b)/(a^3*d) - (a^2 + b^2)*log(e^(-d*x - c) + 1)/(a^3*d) - (a^2 + b^2)*log(e^(-d*x -
c) - 1)/(a^3*d))

Giac [F(-1)]

Timed out. \[ \int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\text {Timed out} \]

[In]

integrate((f*x+e)*coth(d*x+c)^3/(a+b*sinh(d*x+c)),x, algorithm="giac")

[Out]

Timed out

Mupad [F(-1)]

Timed out. \[ \int \frac {(e+f x) \coth ^3(c+d x)}{a+b \sinh (c+d x)} \, dx=\int \frac {{\mathrm {coth}\left (c+d\,x\right )}^3\,\left (e+f\,x\right )}{a+b\,\mathrm {sinh}\left (c+d\,x\right )} \,d x \]

[In]

int((coth(c + d*x)^3*(e + f*x))/(a + b*sinh(c + d*x)),x)

[Out]

int((coth(c + d*x)^3*(e + f*x))/(a + b*sinh(c + d*x)), x)